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Abstract 

This paper presents the design and simulation of full order optimal controller for deregulated power 
sxstem for Automatic Generation Control (AGC), Traditional AGC of two-area system is modified to 
take in to the effect of bilateral contracts on the system dynamics. The DISCO participation matrix 
defines the bilateral contract in a deregulated environment. This paper reviews the main structures, 
configurations, modeling and characteristics of AGC in a deregulated environment and addresses the 
control area concept in restructured power Systems. To validate the effectiveness of full order 
optimal controller, a simulation has been performed using MATLAB and results are presented here. 
The results for LFC and AGC for a deregulated interconnected power systems shows that the 
optimal full order controllers perform better than classical integral order controllers . 

Keywords: Automatic Generation Control, Area Control Error, ACE Participation Factor, 
Bilateral Contracts, Contract Participation Factor. Deregulation, DISCO Participation Matrix, Full 
Order Optimal Controller. Load Frequency Control. 



I. Introduction 

In the electric power system load demand of the consumer always keeps on changing, hence the system 
frequency varies to its nominal value and the tie line power of the interconnected power system changes to its 
scheduled value. AGC is responsible to control the frequency to its nominal value and maintain the tie line 
power to its scheduled value, at the time of load perturbation in the system. In the conventional power system 
the generation, transmission and distribution are owned by a single entity called a Vertically Integrated Utility 
(VIU). In the deregulated environment Vertically Integrated Utilities no longer exist. However, the common 
operational objective of restoring the frequency at its nominal value and tie line power to its schedule value 
remain the same. In the deregulated power system the utilities no longer own generati 
distribution. In this scenario there are three different entities generation companies (GENCOs), I 
companies (TRASCOs), distribution companies (DISCOs). As there are several GENCOs and DISCOs in the 
deregulated environment, a DISCO has the freedom to have a contract with any GENCO for transaction of 
power. A DISCO has freedom to contract with any of the GENCOs in their own area or another area. Such 
transactions are called "bilateral transactions" and these contracts are made under the supervision of an impartial 
entity called Independent System Operator (ISO). ISO is also responsible for managing the ancillary services 
like AGC etc. The objective of this paper is to modify the traditional two area AGC system to take into account 
the effect of Bilateral Contracts. The concept of DISCO participation matrix is used that helps in the 
visualization and implementation of Bilateral Contracts. Simulation of the bilateral contracts is done and 
reflected in the two-area block diagram. The full order optimal controller is used for accomplish the job of 
AGC i.e to achieve zero frequency deviation at steady state, and to distribute generation among areas so the 
interconnected tie line power flow match the prescribed schedule and to balance the total generation against the 
total load. 

II. Formulation of Model of AGC for deregulated power system 

Consider a two-area system in which each area has two GENCOs (non reheat thermal turbine) and two DISCOs 
in it. Let GENCOl GENC02, DISCOl and DISC02 be in area 1 and DISC03, and DISC04 be in area 2 as 
shown in Fig. 1 . 
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For LFC or AGC conventional model is used which is just the extension of the traditional Elgerd model [3]. In 
this AGC model, the concept of disco participation matrix (DPM) is included to incorporate the bilateral load 
contracts. DPM is a matrix with the number of rows equal to the number of GENCOs and number of columns 
equal to the number of DISCOs in the system. The DPM shows the participation of a DISCO in a contract with 
GENCO. For the system described in Fig 1, the DPM is given as 

CP/,1 <?>/l2 <%fn Wfu 

cpfn cpfn cpfx c Pfn 
cpfii <%fn c Pfn c Pfu 
cpf 4l Cjp/ 42 cpf a cgf M 



DPM = 



Where Cpf refers to "contract participation factor." cpfij 



_ Demand of DISCOj from GENCOi 
Total Demand of DISCOj 



Thus ijth entry corresponds to the fraction of the total load power contracted by DISCO j from GENCO i. The 
sum of all the entries of particular column of DPM is unity. 

Whenever the load demanded by a DISCO changes, it is reflected as a local load in the area to which this 
DISCO belongs. As there are many GENCOs in each area, ACE signal has to be distributed among them in 
proportion to their participation in the AGC. "ACE (Area Control Error) participation factor (apf)" are the 
coefficient factors which distributes the ACE among GENCOs. If there are m no of GENCOs then Y^Li a Pfi — 
1 . 

In deregulated scenario a DISCO demands a particular GENCO or GENCOs for load power. These demands 
must be reflected in the dynamics of the system. Turbine and governor units must respond to this power 
demand. Thus, as a particular set of GENCOs are supposed to follow the load demanded by a DISCO, 
information signals must flow from a DISCO to a particular GENCO specifying corresponding demands. Here, 
we introduce the information signals which were absent in the conventional scenario. The demands are 
specified by (elements of DPM) and the pu MW load of a DISCO. These signals carry information as to which 
GENCO has to follow a load demanded by which DISCO. 

The block diagram for two area AGC in a deregulated power system is shown in Fig 2. In this model the 
schedule value of steady state tie line power is given as 

A Ptiei-2,scheduied =(demand of DISCOs in area 2 from GENCOs in areal)-(demand of DISCOs in area 1 from 
GENCOs in area2) 

At any given time, the tie line power error, APtiel-2,error is defined as- 

AP tie l-2,error = AP tiel _ 2 , actuiU - APfei^scheduled 

This error signal is used to generate the respective ACE signals as in the traditional scenario 

ACE , = B ! A/j + &P nel _ 2 errl , r 
ACE 2 = B 2 Af 2 + AP,, r2 _ lem , r 
Where 

P , 

AP = - -^-AP { 

P r2 

P rl and P r2 are the rated powers of area 1 and area 2, respectively. Therefore 

ACE 2 = B 2 Af 2 + a n AP rie ,_ 2 emr 
Where al2 = - — 

Pr2 

In the block diagram shown in figure 2, AP L i iLO c and AP L2 ,loc represents the local loads in area 1 and area 2 
respectively. 

AP L1 . AP L2 . APu and AP L4 represents the contracted load of DISCOl, DISC02, DISC03 and DISC04 
respectively. 
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III. Design of full order optimal Controller 

The theory of optimal control is concerned with operating a dynamic system at minimum cost. The 
case where the system dynamics are described by a set of linear differential equations and the cost is described 
by a quadratic functional is called the LQ problem. The optimal control problem for a linear multivariable 
system with the quadratic criterion function is one of the most common problems in linear system theory, it is 
defined below: 

x=Ax + Bu (1) 

Given the completely controllable plant, where x is the nX 1 state vector, u is the pX 1 input vector. A is the n 
x n order of real constant matrix and B is the nx p real constant matrix. Desired steady -state is the null state 
x=0 

The control law 

u = -Kx (2) 

Where K is pxn real constant unconstrained gain matrix, that minimizes the quadratic performance index . 
The design of a state feedback optimal controller is to determine the feedback matrix 'K' in such a way that a 
certain Performance Index (PI) is minimized while transferring the system from an initial arbitrary state 
x(0) ^ 0 to origin in infinite time i.e., x(oo) = 0 .Generally the PI is chosen in quadratic form as: 

I'l- j/ 0 "(x T Qxfu' R u)dt (3) 

where, 'Q' is a real, symmetric and positive semi-definite matrix called as 'state weighting matrix' and 'R' is a 
real, symmetric and positive definite matrix called as 'control weighting matrix'. 

The matrices A, B, Q & R are known. The optimal control is given by u = - Ajc,'K' is the feedback gain matrix 
given by; 

K = R" 1 B T S (4) 

where, 'S' is a real, symmetric and positive definite matrix which is the unique solution of matrix Riccati 



Equation: 

A T S + SA-SBR- 1 B T S + Q= 0 (5) 

The closed loop system equation is; 

x=Ax + B{-Kx) = (A — BK)x = A c x (6) 



The matrix A — B K — A c is the closed loop system matrix. The stability of closed loop system can be tested 
by finding eigen values of A c . 

IV. State Space Modeling AGC System in Deregulated Environment 

The two area AGC system considered has two individual area connected with a tie line. The deviation in each 
area frequency is determined by considering the dynamics of governor, turbines, generators and load 
represented in that area. The state space model of representation of AGC model is given by 

x = Ax + Bu + tp + Pq (7) 

This model of AGC is shown in Fig. 2 .Where x is state vector, u is control vector p is disturbance vector and q 
is the bilateral contract vector. A, B, X and P are the constant matrix associated with state control, disturbance 
and bilateral contract vector respectively. 

In our system we can identify total 13 states. All these vectors and matrix are given by - 

The State Vector 'x' (13x1), x=[x t x t x 3 x t x i x 6 x 1 x s x 9 x 10 x u x n x n f where 

* = [/, /j AP CV| AP CV , AP GV} AP CV) AP ul A P M , AP (/3 AP M4 j ACE t dt J ACE , dt AP fel _J ; 
(8) 

Control Vector 'u' (2x1) u = [«, uj ; (9) 

Disturbance vector 'p' (2x1) p = [ P <" ; (10) 

Bilateral Contract Vector 'q' (4x1) q = ^- AP ^ AP " A/> " AP "^ (11) 

www.ijceronline.com Open Access Journal Page 30 



Design of Full Order Optimal Controller... 




Fig. 2: Two area AGC model in deregulated power system 
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CP.f 

cpf 41 



cpf i2 
cpf 42 



T 

cpf 2J 



CPf M 

r 

<-'/>/' ; 4 

T 2 

CPf i4 



- (CPf 3, + Cp/ 4 , ) 

- a a (cpf 3 , + cpf 41 ) 



(15) 

V. Design of full Optimal controller for AGC in Deregulated Environment 

The design of a state feedback optimal controller is to determine the feedback matrix 'K' in such a way that a 
certain Performance Index (PI) is minimized while transferring the system from an initial arbitrary state 
x(0) ^ 0 to origin in infinite time i.e., x(oo) = 0. 
Generally the PI is chosen in quadratic form as given by equation (3) 
PI = ijT(x T Qx + u T Ru)dt 

where, 'Q' is a real, symmetric and positive semi-definite matrix called as 'state weighting matrix' and 'R' is a 
real, symmetric and positive definite matrix called as 'control weighting matrix'. The matrices Q and R are 
determined on the basis of following system requirements. 

1) The excursions (deviations) of ACEs about steady values are minimized. In this model, these excursions are; 

ACEi = BjAfi + AP tie !_ 2 = + x 13 and (16) 

ACE 2 = B 2 Af 2 + AP tie 2 _! = B 2 x 2 + a 12 x 13 (17) 

2) The excursions of / ACE dt about steady values are minimized. In this model, these excursions are x n & \ l2 

3) The excursions of control inputs Uj and u 2 about steady values are minimized. 
Under these considerations, the PI takes a form; 

PI = I /"[( B lXl + x 13 ) 2 + ( B 2 x 2 + a 12 x 13 ) 2 + ( Xll ) 2 + (x 12 ) 2 + ( Ul) 2 + (u 2 ) 2 ]dt (18) 
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VI. Simulation Result 

A. Case 1: 

Consider a case where all the DISCOs contract with the GENCOs for power as per the following DPM: 

T .3 .5 .1 .35] 

j .2 .5 .25 .3 
DPM = I 

I .1 0 .4 .1 | 

|_ .4 0 .25 .25 J 
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It is assumed that each DISCO demand s 0. 1 pu MW power from GENCO as defined by cpfs in DPM matrix 
and each GENCO participates in AGC as defined by following apfs: apfi = 0.75, apf 2 = 0.25, apf 3 = 0.5 and 
apf 4 = 0.5. 

The system in figure 3 simulated using the above data and the system parameters given in Appendix -I . The 
result of the simulation has shown in Figure 4. 

The off diagonal elements of DPM corresponds to the contract of a DISCO in one area with a GENCO in 
another area. 



The schedule power on the tie line in the direction from areal to area2 is - 

kPtiel-2,scheduled = i^LS (cp/l 3 + C P/ 23 ) + AP L4 (cp/ 14 + Cp/ 24 )] 

- [AP L1 (cp/ 31 + cp/ 41 ) + AP L2 (cp/ 32 + cp/ 42 )] 

= [0.1(0.1+0.25)+0.1(0.3+0.35)] -[0.1(0.1+0.4)+0.1(0+0)] 

= 0.05 pu MW 

The desired generation of a GENCO in pu MW can be expressed in terms of Contract Participation Factor 
(cpfs) and the total demand of DISCOs as 



AP]v 



iAP Li 



Where AP Lj is the total demand of DISCO j and cpf s are given by DPM. 
At steady state power generated by GENCOs - 

AP M1 = 0.3(0.1) + 0.5(0.1)+ 0.1(0.1)+ 0.35(0.1) = 0.125 pu MW 

AP M2 = 0 . 125 pu MW; AP M3 = 0.06 pu MW; AP M4 = 0.09 pu MW 
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Fig. 4a: Frequency Deviations (Hz) for area 1 & area 2, easel 



Fig. 4a shows the dynamic responses of frequency deviations in two areas (i.e., Afl and Af2).The frequency 
deviation in each area goes to zero in the steady state. 

The schedule power on the tie line in the direction from areal to area2 is .05 pu MW. Fig. 4b shows the actual 
Tie line power flow between two area at steady state is also .05 pu MW. So the deviation in tie line power at 
steady state become zero . 

Fig. 4c & 4d show the power generated by GENCOl, GENC02, GENC03 & GENC04 in steady state. The 
results are matching from our calculations. 
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Fig.4b: Actual Tie line power ( pu MW), case 1 




time (sec) 

Fig.4c: Generated Power by GENCOl & GENCQ2 (pu MW) case 1 
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Fig.4d: Generated Power by GENC03 & GENC04 (pu MW), case 1 
B. Case 2: Contract Violation 

It may be happen that a DISCO violate a contract by demanding more power than that specified in the contract. 
This excess power is not contracted out to any GENCO. This uncontracted power must be supplied by the 
Gencos in the same area as the DISCO. It must be reflected as a local load of the area but not as the contracted 
demand. Now consider DPM as below 

Now in this case all the Discos contract with the Gencos for power as per the following DPM- 



.5 .25 0 


.3 1 


.2 .25 0 


» 


0 .25 1 


.7 1 


.3 .25 0 


1 

° 1 



It is assumed that each DISCO demand s 0. 1 pu MW power from GENCO as defined by cpfs in DPM matrix 
and DISCOl demand 0.1 pu MW of excess power. ACE participation factors are apfl=0.75, apf2=l-apfl= 
0.25, apD= 0.5, apf4=l-apf3=0.5. 

The total local load in area I = Load of DISCOl + Load of DISC02 

= (0.1 +0.1) + 0.1= 0.3 puMW 
Similarly, the total local load in area II = Load of DISC03 + Load of DISC04 

= (0.1 + 0.1) = 0.3 pu MW (no un contracted load) 
Schedule tie line power flow is 0.05 from area 2 to area 1. 

The system in figure 3 simulated again ,using the above data and the system parameters given in Appendix -I . 
The result of the simulation has shown in figure 5. 

Fig. 5a: shows the dynamic responses of frequency deviations in two areas (i.e., Afl and Af2).The frequency 
deviation in each area goes to zero in the steady state. 

Fig. 5b: shows the tie line power deviation between two area. The actual tie line power .05 pu MW which flows 
from area 2 to area l.Scheded tie tie line power is also .05, in the steady state deviation goes to zero. 

Fig 5c & 5d show the power generated by GENCOl, GENC02,GENC03 & GENC04 in the steady state. The 
generation of GENC03 & GENC04 are not affected by the excess load of DISCOl. The un contracted load of 
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DISC01 is reflected in generation of GENCOl & GENC02. The ACE participation factor decide the 
distribution of un contracted load in the steady state. Thus this excess load is taken up by the GENCOs in the 
same area as that of the DISCO making the un contracted demand. 
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Fig. 5a: Frequency Deviations (Hz) for area 1 & area 2, case2 




Fig. 5b: Actual Tie line power ( pu MW), case 2 
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Fig.5c: Generated Power by GENCOl & GENCQ2 (pu MW), case 2 



§ 0.25 




15 



15 



Time (sec) 

Fig.5d: Generated Power by GENC03 & GENC04 (pu MW), case 2 
The value of 'k' optimal controller gain matrix, find out with the help of MATLAB program is given below. 

f .312 -.080 .571 .571 -.128 -.128 .108 .108 -.023 -.023 - 1 0 .147 1 

[-.080 .312 -.128 -.128 .571 .571 .023 -.023 .108 .108 0 - 1 - . 147 J 
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Hence the control inputs: 

ui= -.312 x, +.08 x 2 -.571 x 3 -.571 x 4 + .128 x 5 +.128 x 6 -.108 x 7 -.108 x 8 + .023 x 9 +.023 x 10 +x n - .147 x 13 
ui= .080 x, -.312 x 2 +.128 x 3 +.128 x 4 -.571 x 5 -.571x 6 -.023 x 7 +.023 x 8 -.108 x 9 -.108 X!+x 12 + .147 x 13 

The eigen values of 'A' open loop system are 

[0, 0,-0.444 + 3.631, -0.444-3.63i, -0.823 + 3.021, -0.823 -3.021,-0.778, -13.4, -13.4,-2.50, -2.50,- 
12.5, -12.5] 

Two eigen values are zero and remaining have negative real parts indicating that, the system is marginally stable 
before applying the optimal control strategy. 

The eigen values of 'Ac' closed loop system are 

[-13.3996 -13.3804 -0.8468 + 3.7432i -0.8468 - 3.7432i -1.1232 + 3.1450 -1.1232 - 3. 1450i - 
0.8032 + 0.24351 

-0.8032 - 0.2435i -0.4590 -12.5000 -12.5000 -2.5000 -2.5000] 

All eigen values of 'Ac' have negative real parts indicating that the system is asymptotically stable after 
applying optimal control strategy. 

APPENDIX I 



Prl = P,2 


2000 MW 


R, = R,= R, = R 4 


2.4 Hz /pu MW 


K,= K, 


.6558 




120 


T p i = T p2 




B, = B 2 


.429 


T, 2 


.0707 MW / radian 


T T i =T t: =T t1 =Tt4 




T s i = T s2 =T g ., = T,4 


.08 




200 MW 



VII. CONCLUSION 

This work gives an overview of AGC in deregulated environment which acquires a fundamental role to 
enable power exchanges and to provide better conditions for the electricity trading. The important role of AGC 
will continue in restructured electricity markets, but with modifications. Bilateral contracts can exist between 
DISCOs in one control area and GENCOs in other control areas. The use of a "DISCO Participation Matrix" 
facilitates the simulation of bilateral contracts. Models of interconnected power systems in deregulated 
environment have been developed for integral as well as optimal control strategies. The state equations and 
control equations have been successfully obtained and full State feedback optimal controller has designed .The 
models have also been tested for system stability before and after applying closed loop feedback control and it 
has been observed performance is much better in case of full state feedback optimal controller as compare to 
integral controller. 
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